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ABSTRACT 

We present a study of the kinematics and the physical properties of the central region 
of the Hickson Compact Group 31 (HCG 31), focusing on the HCG 31A+C system, 
using integral field spectroscopy data taken with the Gemini-south telescope. The main 
players in the merging event (galaxies A and C) are two dwarf galaxies that in the past 
have already had one close encounter, given the observed tidal tails, and may now be 
in their second approach, possibly about to merge. We present new velocity fields and 
Ha emission, stellar continuum, velocity dispersion, electron density, Ha equivalent- 
width and age maps. Considering the high spatial resolution of the IFU data, we were 
able to measure various components and estimate their physical parameters, spatially 
resolving the different structures in this region. Our main findings are the following: 

(1) we report for the first time the presence of a super stellar cluster next to the 
burst associated to the HCG 31C central blob, related to the high values of velocity 
dispersion observed in this region as well as to the highest value of stellar continuum 
emission. This may suggest that this system is cleaning its environment through strong 
stellar winds that may then trigger a strong star formation event in its neighborhood, 

(2) among other physical parameters, we estimate an L(H<a)^14xl0 41 ergs -1 and 
a SFR/'Tl M 0 yr -1 for the central merging region of HCG 31 A+C. These values 
indicate a high star formation density, suggesting that the system is part of a merging 
object, supporting previous scenarios proposed for this system. 

Key words: galaxies: abundances, galaxies: interactions, galaxies: kinematics and 
dynamics 


Based on observations obtained at the Gemini Obser¬ 
vatory, which is operated by the Association of Universities 
for Research in Astronomy, Inc., under a cooperative agree¬ 
ment with the NSF on behalf of the Gemini partnership: the 
National Science Foundation (United States), the Science 
and Technology Facilities Council (United Kingdom), the 
National Research Council (Canada), CONICYT (Chile), 
the Australian Research Council (Australia), Ministerio da 
Ciencia e Tecnologia (Brazil) and Ministerio de Ciencia, Tec- 
nologfa e Innovacion Productiva (Argentina) - Observing 
runs: GS-2012B-Q-60. 


1 INTRODUCTION 

Local compact groups of galaxies are excellent places where 
to study galaxy transformation. As a consequence of the 
strong gravitational encounters, galaxies in compact groups 
develop tidal tails, bursts of star formation, gas flows, nu¬ 
clear activity and kinematic perturbations. All these phe¬ 
nomena have strong influences on the evolution of galaxies 
in compact groups. For example, simulations suggest the ex¬ 
istence of gas inflows in merging systems. These flows can 
mix the metal content of the interacting galaxies (Rupke et 
al. 2010a, Perez et al. 2011), and can enhance the star forma¬ 
tion of the central parts of the interacting/merging systems, 
increasing their IR luminosities. On the other hand, such 
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interactions can affect the global kinematics of the galaxies 
in the compact group. This effect can be easily seen in the 
2D velocity fields of such systems (Amram et al. 2003, Plana 
et al. 2003). Given that all the phenomena listed above are 
connected, it will be desirable to analyze all of them simul¬ 
taneously, which in practice is not an easy task, given the 
lack of proper data and/or instrumental limitations. A way 
to circumvent some of these limitations is to use IFU (Inte¬ 
gral Field Unit) Spectroscopic data, given that this tech¬ 
nique allows us to observe simultaneously the spectra of 
several regions, where the galaxy configuration is compact 
enough. We have obtained IFU observations of one of the 
most interesting compact groups of galaxies: Hickson Com¬ 
pact Group 31 (HCG 31, Hickson 1982), in order to study 
the evolutionary phase of this group and to search for a link 
between its kinematics and chemical enrichment history. In 
addition, this system is an ideal laboratory to study star for¬ 
mation at different levels (from small to giant star-forming 
regions), therefore, IFU data on HCG 31 will be extremely 
useful to understand the star formation process in interact¬ 
ing/merging galaxies. 

The Hickson Compact Group 31 is an intriguing object 
because of its wide range of indicators of galaxy interaction 
and merging, like multiple tidal tails, irregular morphology, 
complex kinematics, strong starburst spots and possible for¬ 
mation of tidal dwarf galaxies (TDGs) (Rubin et al. 1990, 
Lopez-Sanchez et al. 2004, Mendes de Oliveira et al. 2006, 
Amram et al. 2007). HCG 31 is formed by several mem¬ 
bers: A+C, B, E, F, G, H, Q and R (Mendes de Oliveira 
et al. 2006). Galaxies A and C form the main central part. 
These are late-type, gas-rich, large Magellanic-type irregu¬ 
lar galaxies (Amram et al. 2007), which could be responsible 
for the two optical tidal tails that extend towards the north¬ 
east and southwest of the A+C system. In addition, TDG 
candidates seem to have been formed out of gas-rich mate¬ 
rial stripped from this system (Lopez-Sanchez et al. 2004). 
Two scenarios have been proposed to explain the nature of 
the system A+C: The system A+C is a single entity, i.e. a 
single interacting galaxy (Richer et al. 2003) and a second 
scenario supports the idea that there are two systems in a 
merging process (Rubin et al. 1990). This latter scenario is 
supported by Amram et al. (2007) who found that the dy¬ 
namics of the A+C system has two main different compo¬ 
nents of approximately the same intensities, indicating that 
it is in a pre-merging stage. 

In this paper we present integral field spectroscopic ob¬ 
servations of the central region of HCG 31. We describe the 
Gemini Integral Field Spectroscopic observations, data re¬ 
ductions and analysis in sections §2 and §3 respectively. Our 
results are presented in section §4. The discussion and con¬ 
clusions are presented in sections §5 and §6, respectively. 
We have adopted a distance of 57.6 Mpc assuming a redshift 
z=0.0137 (Hickson et al. 1992) and Ho = 72kms -1 Mpc -1 . 


2 OBSERVATIONS AND DATA REDUCTION 

The central region of HCG 31 was observed with the 
Gemini South telescope during the nights between Jan¬ 
uary 12 - 19, 2013, with Gemini Multi-Object Spectrograph 
(GMOS, Hook et al. 2004) and the Integral Field Unit (IFU, 
Allington-Smith et al. 2002), under the program GS-2012B- 


Q-60 (PI: S. Torres-Flores). In order to have a larger spa¬ 
tial coverage and map the whole A+C system, three dif¬ 
ferent fields were observed with the IFU in two-slit mode. 
The spectra and calibrations were acquired using the grat¬ 
ing R400+r’-hlter centered at 6300 A. In the two-slit mode, 
the IFU covers a field of view of 7” x 5”. Each IFU hexag¬ 
onal element (the two-slit mode uses 1500 elements) has a 
projected scale of 0.2” on the sky. For each field, three expo¬ 
sures of 700 sec each were obtained (no spectral or spatial 
dithering was applied during the observations). All fields 
were observed during dark time (sky background 50%) at 
an airmass between 1.23 and 1.31, with median seeing val¬ 
ues of ~0.5” and 0.6” and were not taken under photometric 
conditions (cloud cover: 70% cirrus). Figure [l] shows an HST 
optical image (ACS, F435W) with the location of the three 
GMOS/IFU fields superimposed onto it (“east”, “center” 
and “west” fields, respectively). In addition, the spectropho- 
tometric standard star LTT 4364 was observed in January 
15, 2013 to flux calibrate the science spectrum. 

The observations were processed using a modified early 
version of the Gemini/GMOS IRAF package to reduce IFU 
spectrsQ A master bias was generated by combined raw bias 
frames. Then all science exposures, comparison lamps, spec¬ 
troscopic twilight flats and Gemini Calibration Unit lamp 
flats (GCAL flats) were overscan-ed/trimmed, and bias sub¬ 
tracted using the master bias. Fibers identification was per¬ 
formed with the task “gfreduce” using high signal-to-noise 
GCAL flats. The resulting file with the fiber identifications 
was used later as a reference to identify the fibers in the 
comparison lamps, science exposures and spectroscopic twi¬ 
light flats. The GCAL flats were processed by removing the 
GCAL+GMOS spectral response and the uneven illumina¬ 
tion introduced by the GCAL unit. The twilight flats were 
used to determine the relative fiber throughput and to cor¬ 
rect for any illumination pattern in the GCAL flat using 
the “gfresponse” task. Wavelength calibration has been 
obtained using high signal-to-noise CuAr comparison lamp 
using the “GSwavelength” routine. The residual values in 
the wavelength solution typically yielded a rms of ~0.lA. 
The science exposure are then flat-fielded, corrected by the 
relative fiber throughput, wavelength calibrated, sky sub¬ 
tracted and extracted. The same reduction steps were ap¬ 
plied to process the standard star LTT 4364. The sensitivity 
function was derived using the “gsstandard” routine and 
the science exposures were then flux calibrated using the 
routine “gscalibrate” . 

The final step was the construction of data cubes for 
further analysis. Each processed 2D science image was trans¬ 
formed to a 3D data cube (x, y, A) using the “gfcube” rou¬ 
tine. The output data cubes (3 per held) were re-sampled to 
square pixels 0.1” in size and corrected for atmospheric dif¬ 
ferential refraction. To produce a single data cube for each 
held, the three cubes were combined using the program “py- 
MOSAIC” kindly provided by Dr. James Turner (2013, pri¬ 
vate communication). The hnal spectra have a resolution of 
~2.4A in FWHM (measured using the sky line at 6300 A) 
and a dispersion of ~0.68A pix -1 , covering a wavelength 
interval of ~5620A-6980A. 


1 http://drforum. gemini.edu/topic/gmos-ifu-data-reduction- 
scripts/ 
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Figure 1 . Hubble Space Telescope (HST) optical image (ACS, F435W) of HCG 31. The image shows the regions A+C, B and part of 
the southeast tidal tail. The red squares show the location of the three IFUs in the region A+C. These regions are denominated, from 
left to right: region east, center and west. Image published by Torres-Flores et al. (2015). 


3 ANALYSIS 

The optical emission of the central region of HCG 31 A+C 
is dominated by two main bursts of star formation which 
coincide with the centers of galaxies HCG 31A and HCG 
31C (Iglesias-Paramo V Vflchez 1997). The ages of these 
bursts are very similar and therefore it is inferred that they 
were triggered by the interaction between galaxies A and 
C (Iglesias-Paramo V Vflchez 1997). Lopez-Sanchez et al. 
(2004) detected some underlying stellar emission associated 
with these bursts. However, it is negligible compared with 
the strong emission coming from the nebular lines. In this 
work we will focus on the flux coming from the ionized gas. 
In Figure [2] we show a typical spectrum observed in the 
central region of HCG 31A+C. This spectrum is dominated 
by strong emission lines, where the line ratios are typical of 
star-forming regions in the BPT plot (Baldwin et al. 1981), 
with almost no emission coming from an underlying stellar 
population. 

In Figure [3] we have included an optical image of the 
HST /ACS (F435W) for the system HCG 31A+C, where 
we have over plotted the three IFU fields in red rectangles 
(5”x7” or ^1355x1897 pc each). In addition, we have over 
plotted the position of the two main bursts of star formation 
detected by Iglesias-Paramo V Vflchez et al. (1997) in yellow 
squares, which are located in the galaxies HCG 31C and 


HCG 31A (boxes have sizes 1.1” x 1.1”, or ~298x298 pc). 
The box size was chosen considering the burst associated 
to the galaxy HCG 31C, since this burst shows a stronger 
and more extended H<a emission. The mean values of the 
different physical parameters of these two main bursts were 
obtained considering this area. 

3.1 Extinction 

Observed flux should be affected by Galactic and internal 
extinctions. We have corrected the observed flux for Galac¬ 
tic extinction by using a color excess of E(B-V)=0.04 mag 
(NED database) and the Fitzpatrick (1999) extinction law. 
The internal extinction clearly depends on the internal prop¬ 
erties of the observed source. Starburst objects are usually 
dusty and, therefore, we expect to have high internal extinc¬ 
tions. We note that the spectral coverage of our observations 
do not allow us to derive the Ha/H/3 ratio to determine the 
internal extinction through the Balmer decrement. There¬ 
fore, in order to correct for this effect, we have used the 
spectra published by Mendes de Oliveira et al. (2006) to ob¬ 
tain an average nebular color excess for the central region of 
HCG 31 A+C. In this case, the nebular color excess was es¬ 
timated by using the observed H<a/H /3 ratio and comparing 
it with the expected value listed in Osterbrock (1989), us¬ 
ing the recipes given in Dominguez et al. (2013). We obtain 
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an average nebular color excess of E(B-V)~0.08 mag. Given 
that HCG 31A and HCG 31C are starburst objects (Gal¬ 
lagher et ah 2010), we have used the Calzetti et al. (2000) 
extinction law to correct our data by internal extinction, 
considering the nebular color excess estimated above (which 
is a factor of ~2 larger than the stellar color excess). We note 
that the extinction correction is not a critical issue for line 
ratios (given that the emission lines display a small wave¬ 
length separation). However, this correction is important for 
the estimation of Ha luminosities and star formation rates. 


3.2 Flux measurements 

The main nebular emission lines available in our data 
correspond to [NIIJA6548A, HaA6563A, [NIIJA6584A, 
[SIIJA6717A and [SII]A673lA. The available spectral res¬ 
olution of FWHM=2.4 A (~114kms 1 , value adopted from 
method (ii), see next section) does not allow us to resolve 
multiple components on each emission line profile. For this 
reason, fluxes have been estimated by fitting a single Gaus¬ 
sian on each observed line, by using fluxer[^] which allows 
us to determine the continuum level around the emission 
lines in a confident way. The stellar continuum definition 
was made in the wavelength range of 6620-6690 A. The out¬ 
put parameters of this procedure are flux, continuum, width, 
center, equivalent width and uncertainties derived from the 
2D maps. Fluxes were compared with the values estimated 
by fitting a single Gaussian with the mpfitJ^] code and with 
the task SPLOT in iraf, and the results are fully consistent. 
Given that, we have included the Ha emission and stellar 
continuum emission maps in the top and bottom panels in 
Figure [4] respectively. 


3.3 Velocity measurements 

The most intense emission line in our data set is Ha. This 
emission line is usually used to analyze the kinematics of the 
warm gas in galaxies (e.g Amram et al. 2007). We used it to 
derive the velocity held of the central region of HCG 31. As 
described above, the Gaussian fitting process gives us the 
center of each emission line. In this case, the estimated cen¬ 
ter was converted into velocity, yielding a velocity held. In a 
similar way, we used the width of the prohles to produce a 
velocity dispersion map. At hrst order, if we suppose that the 
different emission line prohles can be htted by Gaussians, the 
velocity dispersion a corrected for instrumental and thermal 
broadening could be written as: a = \/— of nst — cr 2 3 ^, 
where a Q bs is the measured velocity dispersion on the ob¬ 
served prohle. The instrumental width ainst was estimated 
in two different ways which resulted in consistent values: (i) 
from the average of the FWHM of the lines from the CuAr 
calibration lamp close to Ha, obtaining a spectral resolu¬ 
tion of ~2.3d=0.1 A, corresponding to ~110kms -1 , leading 
to a resolution R^2850±125 (in Ha), and (ii) following the 
method published by Sanchez et al. (2012), who estimate the 
spectral resolution htting a Gaussian in the most intense sky 


2 Interactive routine in IDL written by Christof Iserlohe. 
http://www.ciserlohe.de/huxer/huxer.html 

3 Fitting program package in IDL written by Craig B. Markwardt, 
NASA/GSFC Code 662, Greenbelt, MD 20770. 


lines in one frame. In our case, we have obtained an aver¬ 
age of ~2.4±0.1 A, corresponding to ~114kms -1 , leading to 
a resolution R~2730±115. Finally, we adopt the value ob¬ 
tained from the second method. From this method we have 
estimated a instrumental width of ainst ~46±10kms -1 , 
where the uncertainty is represented by the standard de¬ 
viation following Sanchez et al. (2012) method. We found a 
change in the value of the instrumental width across each 
observed fields. On average, this variation is of the order of 
~10kms _1 . However, we note that the largest changes in 
this parameters occur at the edge of the fields, while in the 
centre of the fields the variation in the instrumental width 
is close to ~6kms -1 . For the thermal width correction we 
have assumed a <7^=9 km s -1 , that corresponds to a elec¬ 
tron temperature of 10000 K (which is consistent with the 
temperatures estimated by Lopez-Sanchez et al. 2004). The 
uncertainty associated with the corrected velocity dispersion 
was estimated by propagating in quadrature the uncertain¬ 
ties of ai ns t and a Q b s . On top and bottom panels in Figure 
[5] we present the velocity and velocity dispersion maps, re¬ 
spectively. 


3.4 Star formation rates 


There are several tracers of star formation. For example, 
the UV radiation emerging from young hot stars can be 
used to estimate the star formation rate of galaxies. How¬ 
ever, this emission considers stellar populations with ages 
up to ~200 Myrs. Also, this emission is strongly affected by 
dust attenuation. To mitigate these problems, several au¬ 
thors use the UV and IR emission (e.g. Iglesias-Paramo et 
al. 2006 and more recently Eufrasio et al. 2014), in order 
to take into account the emission of dust that is heated by 
the strong UV radiation held. Another way to estimate star 
formation rates is to use the Ha emission, which emerges 
from gas ionized by massive stars. Given the available data, 
the star formation rates in HCG 31 have been estimated by 
using its Ha luminosity and the recipe given in Kennicutt 
(1998), which assumes a continuous star formation process, 
where SFRij a =7.94x 10 -42 xL(H a ). However, we note that 
the SFRs derived from the Ha emission can be underes¬ 
timated given that this emission is affected by extinction, 
which corresponds to one of the most important sources of 
uncertainties in the determinations of SFR when using Ha 
luminosities (Kennicutt 1998). In the case of the current 
analysis, we have derived an average extinction for the re¬ 
gion under study, which comes from the Balmer decrement 
of archival longslit observations. Of course, the extinction 
can vary along the extension of HCG 31 A+C, given the 
strong burst of star formation located in this region. There¬ 
fore, we caution the reader that our SFRs estimations can 
be affected by spatial variations in the extinction, therefore, 


our SFRs can be overestimated (discussed in £5.1). 


3.5 Electron densities 

Our spectroscopic data cover the [SII] lines, [SIIJA6717A 
and [SII]A673lA. We have calculated the electron density 
(N e ) from the observed [SII] A6717 A/ [SII] A6731 A ratio us¬ 
ing the task temden, which belongs to the iraf stsdas neb¬ 
ular package. The task temden implement the five-level 
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Figure 2. Top spectrum: An example of the typical spectrum observed in the central region of HCG 31. The main emission lines are 
labeled on this image. The position of the region which has this spectrum is shown in Fig. [l] (5” x7” centered on RA: 05:01:37.77, DEC: 
-04:15:30.24). Bottom spectrum: Zoom in the continuum emission of the top spectrum. 


atomic model fivel (De Robertis, Dufour & Hunt 1987). 
We are not able to estimate the electron temperatures due 
to the spectral coverage (~5620 A-6980 A). So, in this case 
we have assumed an electronic temperature of T e =10000K. 
This assumption is consistent with the electronic tempera¬ 
ture found by Lopez-Sanchez et al. (2004) for the burst as¬ 
sociated with HCG 31C (T e =10800±300 K, Table 3 in that 
paper). The uncertainties in the estimates are represented by 
the standard deviation. The electron density map obtained 
using this method is shown on the top panel of Figure [6] 

3.6 Her Equivalent Width and age determination 

We estimated the ages of the bursts by comparing the 
measured Ha equivalent width (EW(Ha)) with Starburst99 
(Leitherer et al. 1999) modelled values. Considering the 
fact that the oxygen abundances in the regions of inter¬ 
est range from 12+log(0/H)~8.2 and 12+log(0/H)~8.5 
(Lopez-Sanchez et al. 2004, Torres-Flores et al. 2015), the 
Starburst99 models were generated for a single stellar pop¬ 
ulation with metallicities of Z=0.004 and Z=0.00^] which 
correspond to Z=O.2-O.4Z 0 , and Salpeter Initial Mass Func¬ 
tion (IMF) with masses ranging between 1-100 M 0 . The 
map obtained for the estimated values of EW(Ha) is shown 
in the bottom panel of Figure [6] and the age maps are shown 
in the top and bottom panels of Figure [7] 

4 RESULTS 

In Figures mum and [7| we show different maps that we 
derived from the IFU observations of HCG 31A+C. 

4 We used the conversion Z=Z 0 x iq( 1 ° 9 ( 0 / H)-12+3.1) anc [ as _ 
sumed log(O/H) 0 =-3.1O (Fiorentino et al. 2012), and a solar 
metallicity of Z 0 =0.018 (Allende-Prieto et al. 2001). 


4.1 Ha emission and star formation rates. 

In the top panel of Figure [4] we show the Ha emission flux 
map of the central region of HCG 31A+C. In the figure, 
the black contours represent the stellar continuum emission, 
which was derived from the IFU data and trace the stellar 
emission in this region (see bottom panel of Figure |4j). In¬ 
specting this figure we detect four strong Ha sources, where 
the northern objects are linked with the galaxy HCG 31C, 
while the southern sources belong to the member HCG 31A. 
We also detect diffuse Ha emission across the three observed 
fields. 

The bottom panel of Figure [4] shows the stellar con¬ 
tinuum map. In the case of HCG 31C, there is an offset 
between the peak of the continuum and the Ha emission. 
The continuum peaks 1 arcsec (~270 pc) to the north¬ 
west of the maximum in the Ha emission. This feature 
has not been previously identified in the central region 
of HCG 31A+C. This can bring some new insights re¬ 
garding the star formation processes in this merging sys¬ 
tem. A possible scenario for this phenomenon is discussed 
in the next sections. In general our Ha map is in agree¬ 
ment with the Ha image published by Vflchez & Iglesias- 
Paramo (1998). Using our data set we have derived the 
Ha luminosities of the two main bursts associated with 
HCG 31C and HCG 31 A, which display luminosities of 
Lff a = 1.89x 10 41 erg s -1 and L# a =1.73xl0 41 erg s -1 , re¬ 
spectively. Clearly these values exceed the luminosity of the 
well known Giant HII regions (like 30 Doradus, which has 
a luminosity of Li/ a ~3xl0 39 ergs -1 , Lopez et al. 2014). 
The burst detected in HCG 31A and HCG 31C have been 
studied in the past by Johnson &; Conti (2000), who found 
that the galaxy component AC shows signs of star formation 
over the past ~10 Myr, and a peak in the EW(Ha) distri¬ 
bution corresponding to star formation about ~5 Myrs ago. 
In the case of the star-forming complexes associated with 
HCG 31C, Lopez-Sanchez et al. (2004) detected Wolf-Rayet 
features, which is consistent with a young stellar population. 
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Figure 3. Hubble Space Telescope (HST) optical image (ACS, F435W) of the system HCG 31A+C. From left to right: region east, 
center and west, with a FOV of 5” x7” each. Red boxes indicate the position of the observed IFUs also shown in Figure^ Yellow boxes 
show the positions of the two main bursts. The northern burst associated to the galaxy HCG 31C and the southern burst to galaxy HCG 
31A. Each box has an area of 1.1” x 1.1” (~298x298 pc). 


These authors also suggest an age of 5 Myrs for the complex 
associated with HCG 31C, and an age of 7 Myrs for member 
HCG 31A. However, the current analysis allows us to study 
in detail the physical properties of these bursts, given that 
we have spectroscopic information for each spaxel. 

Following Kennicutt (1998), the above H a luminosities 
were used to estimate the SFRs in the central region of HCG 
31A+C. The right color bar in the top panel of Figure [4] in¬ 
dicates the SFRs in the HCG 31A+C complex. The figure 
shows a peak in the SFR located in the HCG 31A complex, 
which seems to be more compact than the region located to 
the north. We estimate a SFR=1.49 M®yr _1 and SFR=1.37 
M©yr _1 for the northern and southern star-forming regions 
of HCG 31C and HCG 31A, respectively. Then, for the whole 
complex, we estimate a SFR^2.86 M©yr -1 , which is consis¬ 
tent with the value derived by Lopez-Sanchez et al. (2004) 
for the region HCG 31A+C (see their Table 9). When we de¬ 
rived the SFR for each observed field (“east”, “center” and 
“west”), we found that fields “east” and “west” display a 
modest SFR, when compared with the SFR derived for the 
central field (0.69, 9.16 and 1.21 M©yr _1 , respectively). All 
these values are listed in Table 1 including the uncertainties 
associated to these measurements. Integrating the SFR for 
all the observed members of HCG 31A+C, we determine a 
SFR=11.06 M©yr _1 . Gallagher et al. (2010) used UV and 
IR information to derive the SFR of all members of HCG 
31, finding a value of ~10.6 M 0 yr _1 . They note that the 
main contribution in the SFR comes from members A+C 
and E. Some possible reasons for the discrepancies between 
both determinations are mentioned in the §5.1 [ 


4.2 Kinematics: Velocity fields and velocity 
dispersion maps 

In the top panel of Figure [5] we show the velocity field of 
the central region of HCG 31. The uncertainty associated 
to the radial velocities is ~10kms -1 , which was estimated 
adding in quadrature the error in the wavelength calibra¬ 
tion and the FWHM. We found that the main star-forming 
bursts in HCG 31A and HCG 31C display similar radial ve¬ 
locities, which reach values of ^4000 km s -1 . Although we 
have a small field of view in our observations, we observe 
a perturbed velocity field. We can not identify a kinematic 
major axis. This fact suggests that non-circular motions are 
present in the central region of HCG 31, which is expected, 
given the strong signatures of interactions. In general, the 
velocity field presented in Figure [5] is in agreement with the 
kinematic analysis presented by Amram et al. (2007). 

The velocity dispersion map for HCG 31 is presented in 
the bottom panel of Figure [5] The uncertainty associated to 
the velocity dispersion is of 13kms _1 . The two main bursts 
of star formation, belonging to HCG 31A and HCG 31C, 
display velocity dispersions of ~55kms -1 . In this figure, 
the most interesting feature corresponds to a region located 
in HCG 31C, which displays quite high velocity dispersions 
~85kms -1 , which is identified by a blue box. The velocity 
dispersion in this region differs by at least la (13kms _1 ) 
from that of its neighborhood, as calculated by considering 
the largest variation across the field, in an area of 5” x 7” 
(~10kms -1 ). The variation across an area of 2” x 5”, in 
the central region where the feature is located, is smaller 
(~6kms -1 ). Considering this, we can suggest the existence 
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Figure 4. Top panel: Ha map for the central region of HCG 31A+C. In the same figure we have included a color bar indicating the 
estimated SFRs. Bottom panel: Stellar continuum emission for the central region of HCG 31A+C. Both emission maps are in units of 
erg s —1 cm -2 and the black contours represent the stellar continuum emission. 


of a peak in the velocity dispersion located in this area. 
This source lies at the same spatial position as the peak 
in the continuum contours (see bottom panel of Figure [4]). 
Given the medium spectral resolution of our observations, 
this broadening can emerge as the sum of, at least, two com¬ 
ponents of the Ha emission line. These components may be 
associated with ionized gas that is being blown out by a star 
cluster, which can be identified with the continuum informa¬ 
tion. We can now speculate that the offset between the Ha 
and continuum emission at this location (star-forming burst 
in HCG 31C) is produced by a stellar cluster, which is clean¬ 
ing its environment. This process may be triggering star for¬ 
mation in its neighborhood, and that is the reason why we 
observe a strong Ha emitting source close to the continuum 


peak. This kind of phenomenon has been studied in local 
giant star-forming regions, such as Nil (Walborn & Parker 
1992, Barba et al. 2003) and in the Hu region NGC 604 in 
M33 (Malz-Apellaniz et al. 2004, Tosaki et al. 2007, Barba 
et al. 2009, Farina et al. 2012, Martfnez-Galarza et al. 2012). 


4.3 Oxygen abundances and electron Densities 

The oxygen abundance map of the central region of HCG 
31A+C has been recently published by our group (Torres- 
Flores et al. 2015). This map was derived by using the N2 
index, which is defined as: 


N2 = /og([Nn]A6584/Ha), 


(i) 
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Table 1 . Luminosities and SFRs of HCG 31A+C 


Source 

Life/ 1 )* 
Xl0 41 erg s -1 

Life/ 2 )* 

X 10 41 erg s -1 

LffcA 3 )* 
Xl0 41 erg s -1 

SFRW** 

Meyr" 1 

SFR,(2)** 

Meyr” 1 

SFR( 3 )** 

Meyr” 1 

Burst HCG 31C a 

1.43 ± 0.13 

1.55 ± 0.14 

1.89 ± 0.18 

1.13 db 0.11 

1.22 db 0.11 

1.49 db 0.14 

Burst HCG 31A a 

1.31 ± 0.12 

1.42 ± 0.13 

1.73 ± 0.16 

1.03 ± 0.10 

1.12 ± 0.10 

1.37 ± 0.13 

HCG 31 East 

0.63 db 0.06 

0.70 db 0.65 

0.87 db 0.08 

0.50 db 0.05 

0.55 db 0.05 

0.69 db 0.06 

HCG 31 Center 

8.77 ± 0.82 

9.50 ± 0.89 

11.60 ± 1.08 

6.93 ± 0.65 

7.51 ± 0.70 

9.16 ± 0.85 

HCG 31 West 

1.16 db 0.11 

1.26 ± 0.17 

1.53 db 0.14 

0.92 db 0.09 

0.99 db 0.09 

1.21 db 0.11 

HCG 31 All 

10.56 ± 0.83 

11.46 ± 1.12 

14.00 ± 1.10 

8.35 ± 0.66 

9.05 ± 0.71 

11.06 db 0.86 


Considering a box of 1.1” X 1.1”, corresponding to 298x298 pc. 

Estimated using the integrated Ha emission flux. 

Using the measured integrated L Ha luminosity as input to Kennicutt’s (1998) formula (2), and assuming continuous star formation. 


( 1 ) 

( 2 ) 


(3) 


Value before extinction correction. 

Value corrected by Galactic extinction using the Fitzpatrick extinction law. 


Value corrected by Galactic extinction using the Fitzpatrick extinction law, and internal extinction using the Calzetti extinction law 
with a nebular color excess of E(B-V)g as =0.08 mag. 


and using the Marino et al. (2013) calibration. Torres-Flores 
et al. (2015) found a smooth metallicity gradient between 
HCG 31A and HCG 31C (1.5x 10 -4 dexpc -1 ). Inspecting 
high-resolution Fabry-Perot data published by Amram et al. 
(2007), they found a gas flow between both galaxies, which 
confirms the scenario of metal mixing produced by gas flows, 
given the smooth gradient between HCG 31A and HCG 31C. 

In the top panel of Figure[6]we show the electron density 
map of HCG 31A+C, where the contours represent the Ha 
emission. Given the low intensity of the [SII]AA6717,6731 A 
emission lines, most of the information comes from the main 
star-forming burst. The burst associated to HCG 31A dis¬ 
plays an average electron density of N e =228=bl20 cm -3 , 
with a peak of N e ~500 cm -3 . In the case of the galaxy HCG 
31C, it displays an average of N e =82±50 cm -3 , reaching a 
peak of N e ~150 cm -3 . Lopez-Sanchez et al. (2004) esti¬ 
mated the electron density for the different components of 
HCG 31. In the case of HCG 31C, these authors found a 
value of N e =270±70 cm -3 (the value for the HCG 31A was 
not estimated). However, we note that our data allow us to 
obtain the electron density for each spaxel in HCG 31A+C. 
Therefore, this fact can affect any comparison with previ¬ 
ous results, which mainly have been obtained from longslit 
spectroscopy. 


4.4 Ha Equivalent Width and age estimation 

Given that the Ha emission traces recent star formation, 
the EW(Ha) can be used to estimate the ages of these 
events (e.g., Leitherer & Heckman 1995). The Ha Equiv¬ 
alent Width is important since it gives us an estimation 
of the ratio between the ionizing photons coming from 
massive stars and the continuum photons from the stel¬ 
lar population and gas (Cedres, Cepa &; Tomita 2005). 
The bottom panel of Figure [6] shows EW(Ha) values for 
the central region of HCG 31A+C, where the contours 
represent the Ha emission. The EW(Ha) varies between 
18 A<EW(Ha)<501 A. For the two main bursts we have 
estimated a mean value of EW(Ha)= 338±50A for burst 
in HCG 31A and of EW(Ha)= 238=b90 A for burst in HCG 
31C, whose values are similar to those measured in the star 
forming complexes in the “west” region. 


From the EW(Ha) values we have estimated the age of 
the star-forming objects located in the central region of the 
system HCG 31A+C. These maps are shown in Figure [7| In 
the top an bottom panels the ages were estimated by assum¬ 
ing a metallicity of Z=0.004 and Z= 0.008, respectively. In 
both maps the contours represent the Ha emission. Inspect¬ 
ing the top panel (Z=0.004), we can observe ages ranging 
from 5 to 11 (±2) Myrs, through the “east”, “center” and 
“west” regions. In the bursts associated to the galaxies HCG 
31C and HCG 31A, we have estimated a mean age of 6±2 
and 5±1 Myrs, respectively. Considering the uncertainties, 
we can suggest that both bursts have similar ages. The bot¬ 
tom panel in Figure [7] (Z=0.008) shows similar ages with 
small variations, but following the same trend with ages be¬ 
tween 5-10 (±2) Myrs. 

We have also estimated the age at the same location 
of the stellar continuum emission peak (see bottom panel 
Figure [4j) where we obtain an age of 9 Myrs, while the Ha 
emission peaks in the burst associated to HCG 31C display 
an age of 5 Myrs (see top panel Figure [4]). This result is in 
agreement with the scenario we proposed in §4.2| where the 
stellar continuum emission is coming from a stellar cluster 
(~9 Myrs) which would be cleaning the ISM triggering star 
formation around it (~5 Myrs). The lowest ages are found 
in the galaxy HCG 31C and (even considering the uncer¬ 
tainties) are consistent with Wolf-Rayet stars. This fact is 
consistent with the spectrum of this region published by 
Lopez-Sanchez et al. (2004). 

Although previous works have already dated some star¬ 
forming regions in HCG 31, the IFU data allow estimating 
the ages of different structures in a spatially-resolved way, 
like in the star-forming regions in HCG 31A and HCG 31C, 
and the “east”, “center” and “west” regions. Also, the cur¬ 
rent results show that the region where the super stellar 
cluster candidate lies is older than the burst associated with 
the galaxy HCG 31C. This allows us to add a new piece 
of evidence to the scenario we proposed where we have a 
stellar cluster, which is cleaning its environment triggering 
star formation. On the other hand, other authors have esti¬ 
mated the ages for a more global area in this system (e.g. 
Lopez-Sanchez et al. 2004, Johnson & Conti 2000) which 
range from 5 to 10 Myrs. 
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Figure 5. Top panel: Velocity map of the central region of HCG 31A+C. Bottom panel: Velocity dispersion map of the central region of 
HCG 31A+C, where the blue box overplotted represent the feature with high velocity dispersion values (~70-85 km s _1 ). Black contours 
overplotted on the maps represent the Ha emission. 


4.5 Ionizing mechanism 

In order to analyze the ionizing mechanism in the “cen¬ 
ter” region of HCG 31A+C (where most of the Ho emis¬ 
sion arises) we have included an EW(Ha) versus N2 
(log([Nll]/Ha)) diagnostic diagram, which is presented in 
Figure [§] This kind of diagram was proposed for the first 
time by Cid Fernandes et al. (2010), who used it for a gas 
ionization mechanism classification of their galaxy sample, 
considering the limits proposed by Kewley et al. (2001), 
Kauffmann et al. (2003) and Stasinska et al. (2006). Cid Fer¬ 
nandes et al. (2010) point out the utility of these diagrams, 


which arises from the information we can get from both 
physical parameters (EW and N2). First, they mention that 
the line ratio can give information of the ionized gas phys¬ 
ical conditions, and the EW(Ha) provides an estimation of 
the ionization power of the stellar population. Kauffmann et 
al. (2003) present a classification based on the [Nil]/Ha ra¬ 
tio, proposing an upper limit of -0.32<log([Nn]/Ha) to star¬ 
forming galaxies (green dashed line in Figure [8]). Stasinska 
et al. (2006) present techniques to distinguish star-forming 
galaxies from active galactic nuclei galaxies (AGN). One of 
these is also based on the [Nil]/Ha ratio, where the classifica- 
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Figure 6. Top panel: Electron density map for the central region of HCG 31A+C. Bottom panel: Ha equivalent width map for the 
central region of HCG 31A+C. Black and white contours overplotted on the maps represent the Ha emission. 


tion is given by: -0.4<log([Nll]/Ha) for star-forming galaxies 
(red dashed line in Figure [8]); -0.4<log([Nn]/Ha)<-0.2 for 
AGNs (region between red and blue dashed line in Figure 
|8|; and log([Nn]/Ha)>-0.2 for hybrid galaxies (blue dashed 
line in Figure [8]). In these diagrams, the limits proposed for 
the galaxy classification given by Kauffmann et al. (2003) 
and Stasinska et al. (2006) are approximately vertical. This 
fact allows us to carry out this classification considering the 
[Nil]/Ha ratio (Cid Fernandes et al. 2010). 

In Figure [§] we note that the N2 values for HCG 
31A+C are in the range of -1.2 < log([Nn]/Ha) < -0.5 
(purple points), while the EW(Ha) values are in 


1.3 < log(EW(Ha)) < 2.7. In this diagram we can observe 
that the “center” region is placed over a region where the gas 
in the galaxies is ionized by star formation. This corresponds 
to the limits proposed by Kauffmann et al. (2003) (green 
dashed line) and Stasinska et al. (2006) (red dashed line). 
From these results we conclude that the ionizing mechanism 
in HCG 31A+C is star formation, excluding AGN activity. 

4.6 Intensity versus a plots 

In section [T2] we suggest the presence of a super stellar clus¬ 
ter in HCG 31C, which would be cleaning its enviroment and 
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Figure 7. Top panel: Age map estimated assuming an instantaneous star formation law for the central region of HCG 31A+C for a 
metallicity of Z=0.004. Bottom panel: similar to Top Panel, but for Z=0.008. Both cases for a Salpeter initial mass function (IMF) with 
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triggering star formation in its vicinity. This suggestion was 
supported by the measured velocity dispersion at that loca¬ 
tion. This kind of phenomena has been reported for other 
systems, but for the first time, we suggest these events are 
taking place in the main central region of HCG 31A+C, in 
the interface area between the galaxies HCG 31A and HCG 
31C. In order to investigate the presence of this (or maybe 
also other) structure and the plausibility of a possible wind- 
triggered star-forming scenario we have divided the region 
“center” in four zones (considering the major Ha emission 
regions) and we present “cr” versus “flux” diagrams for these 


zones. These are represented by blue squares in the four pan¬ 
els of Figure [9j which are enumerated respectively as Zone 
1, Zone 2, Zone 3 and Zone 4. In this figure, from left to 
right, the panels correspond to the region “center” extracted 
from the optical hst image (Fig. [3]), Ha emission map (Fig. 
[4| with the stellar continuum emission represented by the 
black contours, velocity field map (Fig. [5]) and the velocity 
dispersion map (Fig. [5]), where the contours represent the Ha 
emission. In the last panel of Figure [9] we have overplotted 
a purple square which represents the super shell produced 
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Figure 8 . Diagnostic diagram EW(Ha) versus N2. The values 
presented belong to the region “center”. The limits proposed by 
Stasinska et al. (2006) for the galaxy classification are: red dashed 
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for hybrid galaxies with log([Nll]/Ha)>-0.2. Green dashed line 
represent the limit proposed by Kauffmann et al. (2003) to dis¬ 
tinguish between star-forming galaxies and AGNs, corresponding 
to -0.32<log([Nn]/Ha). 


by the suggested super stellar cluster in this region, with an 
area of 0.58” x0.58” (158x158 pc). 

The “cr” versus “flux” diagnostic diagrams are a useful 
tool to identify the emission line profile broadening mecha¬ 
nism (i.e: Munoz-Tunon et. al 1996, Martinez-Delgado et al. 
2007, Moiseev & Lozinskaya 2012). In the top panel in Fig¬ 
ure [lO] we present the “sigma” versus “flux” diagram for the 
four zones in region “center”, while in middle and bottom 
panels we include the “radial velocity” versus “flux” and “cr” 
versus “radial velocity” diagrams, respectively. In the three 
diagrams, the blue, green, black and red symbols correspond 
to the zones 1, 2, 3 and 4, respectively. Grey symbols repre¬ 
sent the points which do not belong to any of these zones. 
Analyzing the “cr” versus “flux” diagram we can see that the 
Zone 4 (red symbols) shows a different spatial distribution 


in comparison with the others regions, reaching the high¬ 
est velocity dispersions of the region “center” (cr ^60-85 km 
s -1 ) and fluxes up to ~3.7xl0 -15 ergs s -1 cm -2 . 

This difference is also observed in the “radial velocity” 
versus “flux” (middle panel) and “cr” versus “radial veloc¬ 
ity” (bottom panel) diagrams. This zone displays a wide 
velocity range between 3985-4015 km s _1 , with predomi¬ 
nance of high velocity values. This fact could indicate the 
existence of several structures in the region “center”, where 
each structure shows a different kinematic behavior. This 
could suggest that the Ha line profile broadening mecha¬ 
nism in Zone 4 is different to the other zones. On the other 
hand, from the position of the gray symbols in the “cr” ver¬ 
sus “flux” diagram (top panel Figure 10), we can say these 
do not belong to any knot or do not form any shell. These 
are part of the envelope that surrounds Hu regions. 


4-6.1 Inspecting Zone 4 

Considering the distinctive behavior of Zone 4 in the three 
different diagnostic diagrams with respect to the others 
(which suggest a different and more complicated broadening 
mechanism), we have subdivided it in four sub-zones called 
as a, /3, 7 and 5 (see Figure [ 9 ]) in order to analyze with more 
details the phenomena in this zone. This division was made 
in an attempt to determine more accurately the number of 
structures present - while we expect at least two, there may 
be more. 

We have included these three diagnostic diagrams for 
the four sub-zones in the left column of Figure mi where 
the sub-zones a, /3, 7 and 6 are represented by blue circles, 
green triangles, black asterisks and red stars, respectively. 
For clarity, the super shell placed in the sub-zone (3 has been 
highlighted with purple squares in the right hand panels of 
Figure [ll] (corresponding to the purple square in the last 
panel in Figure [9|. 

Analyzing the “cr” versus “flux” diagram (top left panel 
Figure mi we note two independent distributions, which 
are indicated with turquoise arrows. The right arrow indi¬ 
cates a structure placed in the sub-zone /3 (green triangles). 
This presents the highest velocity dispersion values, reach- 
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Figure 10. Top panel: cr versus flux diagram. Middle panel: radial velocity versus flux diagram. Bottom panel: cr versus radial velocity 
diagram. The diagrams consider the zones 1, 2, 3 and 4, which are represented by blue, green, black and red symbols, respectively. Grey 
symbols represent the points that do not belong to any of these zones. 


ing 85kms _1 . On the other hand, the left arrow shows a 
structure which seems to be composed by the other three 
sub-zones of Zone 4. This last structure presents velocity 
dispersions in a range of ~55-80kms _1 . This is not consid¬ 
ering the black asterisks of the sub-zone 7 , which present 
values of 80-85 km s _1 , because these points do not belong 
to the burst associated to the galaxy HCG 31C (see fourth 
panel in Figure |9|. In the second panel in Figure [9j we note 
that the sub-zone f3 presents the highest stellar continuum 
emission (black contours) and lowest Ha emission with re¬ 
spect to the other sub-zones, especially a, where most of 
the burst is located. This fact suggests that in sub-zone ft 
there is an older stellar population than the observed in the 


other sub-zones, specially in the Ha burst. This offset be¬ 
tween the peak in the Ha and stellar continuum emission 
may be linked with star formation triggered by the feedback 
of the super star cluster. The sub-zone a has the most in¬ 
tense Ha emission, of up to ~5x 10 -15 ergss -1 cm -2 , and 
with a values lower than 80kms -1 . Therefore, we can not 
discard the existence of more than one single stellar popula¬ 
tion (SSP) placed in Zone 4. In fact, we report the finding of 
at least two stellar populations, where the oldest is located 
in the sub-zone /?, and the other would be in the other sub¬ 
zones, mostly in sub-zone a, which would be triggered by 
the feedback of the stellar cluster. 
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Figure 11. Top panels: cr versus flux diagram. Middle panels: radial velocity versus flux diagram. Bottom panels: cr versus radial velocity 
diagram. The left column diagrams include the sub-zones in Zone 4, a , /3, 7 and S, which are represented in blue circles, green triangles, 
black asterisks and red stars, respectively. The right column diagrams include the points in the area in the fourth panel in Figure [9] 
(purple square) represented by purple squares, which corresponds to the super shell produced by the possible super stellar cluster in 
sub-zone £ (0.58” x0.58”, 158x158 pc). 


4-6.2 Inspecting the candidate super shell 

In order to locate the points that belong to the super shell 
produced by the super stellar cluster candidate located in 
the sub-zone /3, we have highlighted a region represented 
by a purple square in the last panel in Figure [5] which has 
been plotted independently in the three diagnostic diagrams 
described previously (right panels of Figure [TT|) represented 
by purple squares. Inspecting the “cr” versus “flux” diagram 
in the top right panel in this figure, this structure appears 
as a region delimited by velocity dispersion between ~73- 


86kms _1 and fluxes of 1.0-3.3x 10“ 15 ergss -1 cm -2 , result¬ 
ing in the most intense peak in this figure. Considering that 
the difference in the velocity dispersion values (at least 1- 
sigma~13kms -1 ) between the super shell and its neighbor¬ 
hood, we suggest the existence of a peak in the velocity 
dispersion located where the super shell is located. 

The results found in the “cr” versus “flux” diagram for 
the super shell is consistent with the model proposed by 
Munoz-Tunon et al. (1996). This model describes a shell 
where it is possible to detect the maximum cr value and 
the minimum surface brightness in its center. While the 
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distance to the center is growing till the internal edge of 
the shell, the surface brightness increases and a decreases. 
This super shell would correspond to the region delimited 
as the possible super shell (purple squares). Analyzing the 
“flux” versus “radial velocity” diagram in the middle right 
panel in Figure m we can see a vertical pattern in the 
purple squares which represent the super shell in the sub¬ 
zone p. This pattern describes a wide range of radial ve¬ 
locities (3985-4010 km s -1 ) in a small range of intensities 
(1.0-2.2x 10 -15 ergs s -1 cm -2 ), which may indicate a local 
expansion of the gas. Finally, in the “radial velocity” ver¬ 
sus “cr” in the bottom right panel in Figure the purple 
squares distribution does not present a clear pattern. This 
could indicate that the expansion is isotropic, with no de¬ 
fined direction. This kind of diagram and analysis (and pat¬ 
tern) was also performed by Bordalo, Plana & Telles (2009), 
who found similar results for specific areas in an Hll region 
in the galaxy Zw 40. 

The ages estimated for the sub-zones a and f3 are con¬ 
sistent with the scenario proposed above. The estimates (see 
£4.4) give an age of ~9±2Myrs for the super stellar cluster 
in sub-zone /3, and of ~5=b2 Myrs for the Ha emission peak 
in sub-zone a , obtaining a difference in age of ~4 Myrs. 


5 DISCUSSION 

5.1 The main physical properties of HCG 31A+C 

The Hickson Compact Group 31 has been studied by sev¬ 
eral authors, in different wavelengths, over the last years. 
Rubin et al. (1990), Williams et al. (1991), Iglesias-Paramo 

6 Vilchez (1997), Vilchez & Iglesias-Paramo (1998), John¬ 
son & Conti (2000), Verdes-Monetengro et al. (2005), Richer 
et al. (2003), Amram et al. (2004, 2007), Lopez-Sanchez et 
al. (2004), Mendes de Oliveira et al. (2006) and Gallagher 
et al. (2010) studied the different physical processes that 
are taking place in this object. Several of these results sug¬ 
gest that the central region of HCG 31A+C is dominated by 
young and strong bursts of star-format ion, which have ages 
lower than 7 Myrs. Rubin et al. (1990) and Lopez-Sanchez 
et al. (2004) reported the presence of Wolf-Rayet features 
in the spectra of this system, which is consistent with even 
younger ages. Also, Mendes de Oliveira et al. (2006) found 
that the galaxy HCG 31C follows the luminosity-metallicity 
relation (in the K-band) defined by dwarf irregular galaxies. 

In this work we have estimated a luminosity 
in Ha for the region “center” in HCG 31A+C of 
log(L# a )=42.06 ± 0.09 erg s -1 . This value is higher than 
the estimate done by Iglesias-Paramo &; Vilchez (1997), who 
also studied the H a emission in HCG 31, obtaining a lumi¬ 
nosity in Ha of log(Lij a )=41.66 erg s -1 for the system HCG 
31A+C. In order to correct by internal extinction, these au¬ 
thors used a colour excess of E(B-V)=0.15. Lopez-Sanchez 
et al. (2004), using the data published by Iglesias-Paramo 
V Vilchez (1997) estimate an Ha luminosity for this system 
of log(Lif a )=41.78 erg s -1 and of log(L# a )=41.54 erg s -1 
from their own spectroscopic data. The differences among 
these values can arise from the color excess and the extinc¬ 
tion law used. In our case, we estimate an average value for 
the color excess, which was derived from the longslit ob¬ 
servations of the system HCG 31A+C. This correction was 


applied to the three observed FOV. In this sense, our lumi¬ 
nosities can be over/under estimated, given the use of an 
average value for E(B-V). Also, we need to keep in mind 
that the dust distribution in this system is not uniform. In 
fact, Gallagher et al. (2010) observed that the dust emis¬ 
sion is higher in the interface zone between galaxies HCG 
31A and HCG 31C. Therefore, considering we have used the 
slit which cross this region for the color excess estimate, 
this could over estimate the fluxes in the rest of the fields 
observed. However, the main part of this study is focused 
in the interface region, where the slit is located, therefore 
this could work at least in this context. Other facts to be 
considered in the different luminosity values are the flux cal¬ 
ibration, which is not absolute and could result in an over 
estimation of the fluxes, the area over which we have in¬ 
tegrated the luminosities and the observational technique 
(IFU vs longslit vs Ha imaging). 


The velocity field of HCG 31A+C is perturbed. Clearly, 
even if the grand pattern motion displays globally a rota¬ 
tion movement, it does not show the fingerprint of a simple 
disk in rotation, given the complex structure of the system. 
This fact is consistent with the map derived by Amram et 
al. (2007), which covers the entire compact group. On the 
other hand, the spatial resolution of the IFU observations 
allows distinguishing structures in the central region of HCG 
31A+C, which present different velocity dispersion values. 
This fact allowed us to detect a structure with the highest a 
values in the central region of HCG 31A+C, which is located 
in one of the “star-forming complexes” (complex 1) detected 
by Gallagher et al. (2010, their Figure 4). These authors pro¬ 
pose that these complexes could be formed by several stellar 
clusters, so the complex 1 could include the star formation 
burst associated to the galaxy HCG 31C and a stellar clus¬ 
ter, which we suggest would be cleaning its environment and 
triggering star formation in its neighborhood. 


Recently, Krabbe et al. (2014) estimated the electron 
density for a few Hll regions located in interacting galax¬ 
ies. These authors found that Hu regions in interacting 
galaxies display systematically higher electron densities than 
those derived for Hu regions belonging to isolated galax¬ 
ies (taken from the literature). These authors found that 
the mean electron densities of star-forming regions in in¬ 
teracting galaxies are in the range of N e — 24-532 cm -3 , 
while those obtained for isolated galaxies are in the range of 
N e = 40-137 cm -3 . Our estimations of electron density for 
HCG 31A and HCG 31C are in the range displayed by the 
interacting sample of Krabbe et al. (2014), which is consis¬ 
tent with the interaction scenario for HCG 31A+C. 


For the first time, we report the detection of a super 
stellar cluster candidate which could be an evidence of cur¬ 
rent and triggered star formation process in HCG 31A+C. 
In addition, despite several studies have been developed to 
understand the compact group HCG 31, the spatial reso¬ 
lution of the IFU observations allowed us to carry out a 
detailed analysis of different structures in the system HCG 
31A+C. With this information we estimate several physical 
properties in this system and we describe the star formation 
bursts in this region in a more detailed way. 
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Table 2. Main physical properties of bursts in HCG 31A+C 


Source 

12+log(0/H)( 1 ) 

N e ( 2 >* 

cm -3 

EW(Ha)< 3 )* 

A 

Age( 4 )* 

Myr 

Age( 5 )* 

Myr 

Burst in HCG 31C a 
Burst in HCG 31A a 

8.22 ± 0.17 
8.44 ± 0.16 

82 ± 50 
228 ± 120 

238 ± 90 

338 ± 50 

6 ± 2 

5 ± 1 

5 ± 2 

5 ± 1 


Considering a box of 1.1” X 1.1”, corresponding to 298x298 pc. 

Mean value, where the associated uncertainties correspond to the standard desviation. 

^ ^ Estimated by Torres-Flores et al. (2015) using the N2 calibrator published by Marino et al. (2013). 

( 2 ) 

Estimated using the task temden in the IRAF STSDAS nebular package. 

(3) 

Obtained with using the idl routine fluxer written by Christof Iserlohe. 

(4) 

Using the measured EW(Ha) as input to the Starburst99 code, assuming instantaneous star formation with 
a metallicity of Z=0.004 and a Salpeter IMF (Leitherer et al. 1999). 

(5) 

Using the measured EW(Ha) as input to the Starburst99 code, assuming instantaneous star formation with 
a metallicity of Z=0.008 and a Salpeter IMF (Leitherer et al. 1999). 


5.2 HCG 31 as a pre-merger object 

Because of its complex morphology, several authors have 
suggested different scenarios to explain the origin of HCG 
31. As mentioned earlier, Richer et al. (2003), using Fabry- 
Perot and spectroscopic data, suggest that galaxies HCG 
31A and HCG 31C are a single entity. Oppositely , Amram et 
al. (2007) used high resolution Fabry-Perot data to suggest 
that both systems are independent late-type galaxies, in a 
pre-merger phase. This last scenario is also supported by the 
results presented in this paper. 

In this context, the SFRs estimated can give us evi¬ 
dence about the merging scenario supported by Amram et 
al. (2007). Considering the area for the SFR estimation in 
the region “center”, which correspond to ~2.57 kpc 2 , we ob¬ 
tain a star formation density of logE=0.55 M® yr _1 kpc -2 
for a SFR^9.2 M® yr -1 . Although we do not have the infor¬ 
mation about the quantity of atomic and molecular gas for 
HCG 31A+C, this star formation density locate this object 
in the upper region of the logS^ as versus logE sfr plot pre¬ 
sented by Daddi et al. (2010, Figure 2 in that paper). This 
fact suggest that HCG 31A+C is a strong starburst object. 
It will be difficult to reconcile this high star formation den¬ 
sity with a single low mass galaxy, therefore, these results 
support the scenario on which HCG 31 is a merging object 
and not a single entity. 

5.3 Current and triggered star formation inside 
HCG 31A+C 

In §4.2| and §4.6| we suggest the existence of a super stellar 
cluster located ~1” (270 pc) from the star formation burst 
associated to the galaxy HCG 31C. This super stellar clus¬ 
ter candidate presents an intense stellar continuum emission, 
with respect to its vicinity. Besides, this object is spatially 
associated with a region with high velocity dispersions up to 
~85kms -1 . This fact suggests that the super cluster would 
be cleaning its environment through a shell in expansion 
which should be produced by the strong stellar winds and/or 
the energy coming from supernovae events. The strong in¬ 
teractions of these winds and the environment seems to be 
triggering star formation in its neighborhood, since we can 
detect a strong Ha emission 270 pc away. In order to check if 
this scenario is possible, we can estimate the velocity of the 
expanding bubble associated with the super stellar cluster, 


in order to determine if this quantity is in the range of ex¬ 
pected values. Considering that the distance from the stellar 
cluster to the Ha peak is about ~270 pc, and taking into 
account that the age of the complex is about 5Myrs (this 
work and Lopez-Sanchez et al. 2004), we derive an expand¬ 
ing velocity of ~53kms -1 . This value is consistent with the 
expanding velocity of some giant Hu regions located in the 
Large Magellanic Cloud (LMC). For example, Rosado et al. 
(1996) measured an expanding velocity of 45kms _1 in the 
central region of Nil. Clearly, the star-forming complex in 
HCG 31 is much brighter than regions like Nil, even 30 
Dor, however, this estimate does not discard the scenario 
proposed above. In fact, triggered star formation has been 
found in local galaxies, like the LMC, where it is possible to 
resolve individual stars (see review by Elmegreen 2011 on 
this topic). Other region widely studied is NGC 604, in the 
spiral galaxy M33. 

The proposed scenario to explain the presence of the 
super stellar cluster candidate, which is triggering star for¬ 
mation in HCG 31C, is consistent with the one proposed 
by Tosaki et al. (2007) for the giant Hu region NGC 604. 
These authors detected a high CO(J=3-2)/CO(J=l-0) ra¬ 
tio gas with an arclike distribution (called by the authors 
“high-ratio gas arc”) around the stellar cluster of NGC 604. 
These authors proposed a relation between the “ high-ratio 
gas arc ”, the stellar cluster and the “new” stars in the Ha 
shell, suggesting the following scenario: In a first place, there 
would be what they called a “first-generation star forma¬ 
tion” associated to the central stellar cluster. Then, by some 
common phenomenon such as stellar winds and/or super¬ 
nova explosions, there should be compression of the sur¬ 
rounding ISM which would result in a dense gas layer formed 
like an arclike distribution (the “ high-ratio gas arc”). Then, 
new stars would be formed within the dense gas layer, what 
they called a “second-generation star formation”, which is 
triggered by the stellar cluster (the first-generation stars). If 
this proposed scenario is correct, one could then observe this 
second star population in the mentioned dense gas layer in 
two ways, as shells emitting in Ha and compact radio con¬ 
tinuum sources. Tosaki et al. (2007) mentioned that these 
facts suggest that the Hu regions in the “high-ratio gas” arc 
are regions with current star formation, and the regions with 
past star formation corresponds to the central star cluster, 
which has dispersed the ISM around. The scenario proposed 
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by Tosaki et al. (2007) for NGC 604 was supported by other 
authors (e.g. Martmez-Galarza et al. 2012, through the use 
of infrared data). In order to confirm this scenario in the 
central region of HCG 31A+C we would need new infrared 
and CO information. 

In addition, the scenario presented in this work is con¬ 
sistent with the “two-stage starburst” model proposed by 
Walborn & Parker (1992). In this model there is a concen¬ 
trated star-forming burst which triggers a second burst af¬ 
ter 2 Myrs. This second generation of stars is formed in the 
molecular clouds located in the periphery of the primary 
star-forming burst, which would be caused by the energy 
coming from massive stars. We found a difference of 4 Myrs 
between the super star cluster and the star-forming burst in 
HCG 31C, and considering the uncertainties, the scenario 
we suggest is consistent with this model. This star-forming 
mechanism is described by Elmegreen (1998), who called it 
“large scale triggering: shells or rings”. In this mechanism 
the ionized gas expansion triggers shell formation, which 
triggers star formation in its periphery. Elmegreen (1998) 
points out that the shell size depends on the pressure and 
environment density, where in low density environments the 
shell can reach a size of a few hundred parsecs. These values 
are reached by the super shell, which we suggest is triggering 
star formation in HCG 31A+C (~300pc). 


6 CONCLUSIONS 

For the first time, Integral Field Spectroscopy data of the 
merging system HCG 31 is presented, specifically of HCG 
31A+C. With all the data analysis and the different esti¬ 
mates of physical properties for this system, we suggest a 
triggered star formation scenario, where we propose the ex¬ 
istence of a super stellar cluster candidate next to the star 
forming burst associated with HCG 31C. The strong stel¬ 
lar winds associated with this source seem to be triggering 
a star formation event in its neighborhood. Our data set 
supports the scenario that suggests that the central region 
of HCG 31 is a merger between two late-type spiral galax¬ 
ies. This merger process has triggered the star formation 
in the center of HCG 31, namely, HCG 31A+C. The study 
we present in this paper shows the power of IFU data in 
analyzing interactions in merging galaxies. 
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